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Properties of the weakly non-ideal Bose gas are considered without suggestion on C-number 
representation of the creation and annihilation operators with zero momentum. The "density- 
density" correlation function and the one-particle Green function of the degenerated Bose gas are 
calculated on the basis of the self-consistent Hartree-Fock approximation. It is shown that the 
spectrum of the one-particle excitations possesses a gap whose value is connected with the density 
of particles in the "condensate". At the same time, the pole in the "density-density" Green function 
determines the phonon-roton spectrum of excitations which exactly coincides with one discovered 
by Bogolyubov for the collective excitations (the "bogolons"). 

PACS number(s): 05.30.Jp, 03.75.Kk, 03.75.Nt, 05.70.Fh 



I. INTRODUCTION 

It is well known that theoretical description of the superfluidity phenomenon must take into account the interaction 
between helium atoms. Also, the closeness of the transition temperature in superfluid state T\ to the temperature To 
of Bose condensation for the ideal Bose gas and accumulation of the macroscopic quantity of particles in the state 
with the momentum zero ("condensate") [1] are the most important features which permit us to apply the models of 
weakly non-ideal Bose gas to the liquid Hell. Therefore, it seems possible that the phenomenon of superfluidity can 
be described within the framework of the weakly non-ideal Bose gas model. 

It is necessary to stress here that the model of the ideal Bose gas does not satisfy the Landau criterion of superfluidity 
[2,3] and cannot explain the specific behavior of the thermodynamic properties of the superfluid helium [4]. However, 
the application of the standard perturbation theory expansion on the inter-particle interaction to the Bose gas at 
T < To at once faces the problem of the appropriate description of the condensate (see, e.g., [5]). 

Therefore, the development of the microscopic theory of the degenerated Bose system can be conducted in two 
ways: 

A) The formulation of some special suggestions (conditions) for some functions (operators) for the degenerated 
Bose gas, 

B) Using another initial model instead of the ideal Bose gas for developing the perturbation theory. 

Starting with the classical Bogolyubov's papers [6,7], the microscopic theory of the degenerated Bose gas has been 
based on the special suggestion on the C-number representation of the creation and annihilation ao operators of 
the particles with the momentum p = 0. 

The results obtained by Bogolyubov, including the spectrum of the collective excitations (that further be called 
"bogolons" after the name of the author), permit us to give the qualitative explanation of the experimental data 
in superfluid helium and to satisfy the Landau condition of superfluidity. It should be mentioned here that the 
mathematical methods of the quantum field theory were first applied to the study of the non-ideal degenerated Bose 
gas by Belyaev [8], without special suggestion on the C-number representation of the operators and ao. In [8], the 
special diagram technique for the perturbation row at zero temperature, was developed and then generalized for the 
one- particle Green functions. In particular, it was suggested to consider (apart from the usual Green function with 
one incoming and, respectively, one outgoing external lines), the additional Green functions with two incoming and 
two outgoing external lines. However, due to the complexity of its mathematical approach, the method suggested in 
[8] is not widely used. 

Hugenholtz and Pines in [9] reformulated the problem at the beginning, by changing, according to Bogolyubov, 
the operators Oq and ao by C-numbers, and they could use almost automatically the quantum field theory methods 
for investigating the Bose gas with the "condensate". In particular, in [9], it was shown that for the C-number 
representation of the operators aj and ao in the excitation spectrum connected with one-particle Green function, the 
gap cannot exist. Their approach is used presently in the theory of the degenerated Bose gas (see, e.g., [5]). 

However, rigorous proof of the correctness of the C-number representation of the operators a,Q and ao is not possible 
and, due to this circumstance, the correspondence of the initial Hamiltonian to the Hamiltonian which arises after 
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changing the operators Oq and ao on C-numbers, remains undecided [10,11]. Moreover, in [12-14], the attempts were 
made to suggest some canonical transformation of the field operators, which would not be connected with the C- 
number representation of the operators Oq and ao (and therefore alternative to the Bogolyubov's one). 

Importantly, in [13-16], the possibility of the existence of two different spectra simultaneously - the one-particle one 
with a gap and the collective one - corresponding to the Bogolyubov's branch of excitations, has been suggested and 
discussed. 

In addition, in [14,18], most of the known results for the thermodynamic properties of the degenerated weakly 
non-ideal Bose gas have been reproduced by using the "dielectric formalism" without C-number representation of 
the operators Oq and ao. Furthermore, in [19], it was argued that the formal restrictions on the use of the standard 
temperature diagram technique for the temperatures T < T are absent. 

Another essential circumstance, connected with the C-number representation of the operators aj and a , has to be 
mentioned. The C-number representation of these operators is based on the commutation relations and the smallness 
of the parameter 1/Nq, where Ao is the operator of the particle number in the state p = 0. Therefore, this parameter 
makes sense only for the system with fixed number of particles, i.e. in the framework of the canonical ensemble where 
No is equal to the average < Ao >. Although the Bogolyubov's canonical transformation of the initial Hamiltonian 
is realized in the canonical ensemble, the calculations of the averages, usually are executed in the grand canonical 
ensemble. 

In the present paper we are exploring the second way, B, for developing the microscopic theory of the degenerated 
Bose system. We are take into account not only the provision of the Landau superfluidity condition, but also the 
consecutive description of the weakly non-ideal Bose gas for the transition from temperatures T > To to the case of 
T < T including the case of a strong degeneration. As an initial approach, we consider the self-consistent Hartree-Fock 
approximation, which is the best one-particle approximation for the normal systems (see, e.g., [20]). 



II. HARTREE-FOCK APPROXIMATION FOR 5(r) - POTENTIAL AND THE GAP IN THE 

ONE-PARTICLE SPECTRUM 

Let us consider the weakly non-ideal Bose gas which consists of the particles with the mass m and zero spin. The 
interaction between the particles is described by the potential U(r) with the Fourier-component u(q) 

lim u(q) = tt(0) = u(q = 0) > (1) 

q^Q 

Expression for the one-particle distribution function /(p) reads 

/(p)=<a+a p >, n=— 5^/(p), (2) 

v 

where a+ and a p are the creation and annihilation operators for the particles with the momentum hp, n =< N > /V 
is the average density of the particles in the volume V at temperature T, N — J2 P a p a pi an d the angle brackets <...> 
denote the averaging in the grand canonical ensemble with the chemical potential ii. 

Now let us pay attention to the spectral representation which shows that the function /(p) is not negative for the 
arbitrary wave vectors p 

/(P) > (3) 

In the framework of the self-consistent Hartree-Fock approximation for T > To, the distribution function /(p) satisfies 
the relation (see, e.g., [21]) 

m= { exp (m^y 1 y 1 (4) 

where E(p) is the energy of the one-particle excitations 

E(p) = s( P ) + nu(0) + u (l)f(P + <0 ( 5 ) 

e(p) = h 2 p 2 /2m is the energy spectrum of a free particle, the second and the third terms on the right side of © 
correspond to the inter-particle interaction in the Hartree and in the Fock (the so-called exchange interaction, which 
is conditioned by the identity of the particles) approximations, respectively. 
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As it will be clear from further consideration, the important point is connected with the presence in ® of the 
distribution function /(p) determined by (@]), but not the distribution function f ld (p) of the ideal gas with spectrum 
of the free particles, as it is usual in the perturbational Hartree-Fock approximation. Therefore, equations (H,© form 
the closed system of equations for determination of the distribution function and the one-particle excitation spectrum 
for the fixed thermodynamic parameters. In the language of the temperature diagram technique [5,21], it means the 
exact summation of some class of the diagrams in the equation for the one-particle Green function. 

For the particular case when the interaction potential has the form of the (^-potential 

U(r) = u(0)6(r), u(g)=tt(0) (6) 

from (@]) and © we obtain 



e(p) - 7(p) - M* 



-l 



/(p) \ < ( — J 11 — ) i (7) 

//= M -2mx(0), 7 (p) = ^)/(p) (8) 

As is follows from Eqs. (Ill),® before transition to the thermodynamic limit, it is necessary to take into account that 
the system is located in the finite volume V . The importance of this point has been mentioned in [9] in the case of 
the ideal Bose gas condensation. 

Taking into account ®, one can see that Eqs. ©,® make sense only for fi* < 0. As this takes place, the 
contribution of the function 7(p) in ©,® is negligible for the arbitrary values of the vector p. In this case, the 
distribution function f(p) is equivalent to the one for the ideal Bose gas f id (p) with the chemical potential ^ changed 
to /i*. Correspondingly, in the transition from /i* to the chemical potential \i %d , all known results for the ideal Bose 
gas are reproduced up to the temperature T of Bose condensation, including the transition temperature itself. 

However, the situation changes dramatically for the temperature of Bose condensation To, when the accumulation 
of the particles in the state with the momentum p = starts, and also below this temperature. 

By analogy with the case of the ideal Bose gas (see, e.g., [1]) it would seem that one can suppose that 

/i* = 0, or^i = 2rai(0), (9) 

and that the function /(p) can be represented in the form 

/(p) =< N Q > S pfi + / T (p)(l - S p , ), (10) 

where A = a^ao is the operator of the quantity of particles with the momentum equal zero ("condensate"), / T (p) 
is the one-particle distribution function with non-zero momenta (the "overcondensate" states) and 

1 " (S) } ' /T(P) = /S(P) = { 6XP ( T ) " (U) 

It is necessary to mention that the relation /(p) =< Ao > S p q (see Eq. (flOl) for / T (p) = 0) was first suggested in 
[22]. 

However, the representation |(9|) for the chemical potential fi is twice as big as the chemical potential of the ideal 
Bose gas in the case of strong degeneration. Besides, the function /y(p) possesses the singularity at small wave vectors 
P- 

Let us mention now that neglecting the value 7(p) at p = for the temperatures T < T , is not correct for 
calculating the distribution function /(p) in ((7]) 



7 (p)=n u(0), no = ^^ (12) 

Therefore, if representation lflQ|) is true, the Eqs. ([9|). lflT|) are wrong. Taking into account (flO|) . lfl2|) from Eqs. 10,® 
for T < To one directly finds 

A* = (2n-no)u(0), (13) 

f T (p) = jexp^^ -lj (l-<5 p ,o), E*(p)=e(p)+n o u(0), (14) 
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On the basis of (fT3l) - lfl5|) it is easy to establish the correctness of the following relations 

lim nn = n, lira u = nii(O), (16) 

T^O T->0 

A = lim£*(p) =n u(0), (17) 

Hm/ T (p) = |exp^-l| <oo, Kmf T (p)=0. (18) 

Therefore, in the framework of the self-consistent Hartree-Fock approximation (@]),((5]) for the one-particle distribution 
function f T (p), for the temperatures T < Tq we can establish that: 

(A) In the case of strong degeneration, the obtained results (fl6| coincide with the known relations for the weakly 
non-ideal Bose gas [5]; 

(B) In the one-particle excitation spectrum, the gap ifTTj) arises between the "condensate" and the "overcondensate" 
states. The appearance of the gap (fl7f is conditioned by the existence of the "condensate". The spectrum of 
the one-particle excitations, taking into account the "condensate" with the zero energy and the spectrum for the 
"overcondensate" states, satisfy the Landau condition of superfluidity; 

(C) The distribution function for the "overcondensate" states is finite at small values of the wave vectors distinct 
from the distribution function l|lip for the ideal Bose gas. On the basis of the above, it is possible to assert that 
the self-consistent Hartree-Fock approximation for the one-particle distribution function |(4|), © is suitable as an 
initial approximation for constructing the theory of Bose gas that takes into account the interaction at arbitrary 
thermodynamic parameters. It is necessary to mention further that according to 1|14|1 in many applications, the 
condition T — * is equivalent to the condition 

r«A. (19) 



III. COLLECTIVE EXCITATIONS AND THE DIELECTRIC FORMALISM 



Having the above results, let us consider the problem of the collective excitations in the weakly non-ideal Bose gas on 
the basis of the "dielectric formalism" [17,18]. The experimental determination of the collective excitations spectrum 
is interpreted on the basis of the existing data on the well observed maximums [23] in the dynamical structure factor 
S(q, lu) for q ^ 0, 

i r°° 

S(q,uj) =— dtexp{iut) < p g (t)p_,(0) >, (20) 

* J — oo 



Pq (t) = <- q/2 (*) a p+q/2(*)> (21) 

p 

where p q (t) is the Fourier-component of the operator of particle density in the Heisenberg representation. The 
dynamical structure factor ifSOj) is directly connected [24] with the retarded density-density Green function x H (q, z ) 
which is analytical in the upper semi-plane of the complex z (Imz > 0) , 

S(q,o>) = —, f z m W*(g,o; + zO), (22) 

1 — exp(—fku/l ) 



X R (q,z) - -jry j o dtexp{izt) < [p, (*)/>_, (Q)] »= - « Pq | p- q » 2 , (23) 

The equalities (j20"l) . (|22|) have to be taken in the thermodynamic limit: V — > oo, < N >^ oo and < N > /V — > const. 
Eq. l(22|) serves as the basis for calculation of the function S(q,w) for quantum systems by the perturbation methods 
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of the diagram technique [5,21,24]. The retarded Green function x R (Qi z ) (|23| is the analytical continuation of the 
temperature Green function x T (q> i^n) 

X T {q,ittn) = — «p q \ P- q »iO n , (24) 

from the discrete multitude of the points on the imaginary axis ifl n = i2imT to the upper semi-plane of the complex z 
[5,21,24]. For the function x T (g, ifin), there exists the diagram representation which is connected with the extraction 
of the irreducible (in g-channel on one line of interaction u(qj) part of x T (g, i£l n ) - the so-called the polarization 
operator H(q,itt n ) [21]. After analytical continuation of the function x 7 (q,^n), we arrive at the expression 

X « M = «. (25) 
s(q,z) 

Here, the function e(q, z), by analogy with the terminology accepted in the theory of Coulomb systems [24], is called 
dielectric permittivity 

s{q,z) = l-u(q)IL(q,z). (26) 

It should be noted that all relations mentioned above in this section are valid for arbitrary interaction u(q). 

Determination of the appropriate approximation for the polarization operator permits us to find the poles of the 
Green function x R (q, z )- These poles describe the collective excitations in the system, which are the solutions of the 
equation 

e(q,z) = 0. (27) 

The above equation, meanwhile, is well known from the theory of the Coulomb systems [25]. 

In considering the case of the weakly non-ideal Bose gas for calculation of the polarization operator H(q,z), we 
restrict ourselves to the simplest "one-loop" approximation, which in the theory of the Coulomb systems [25] is called 
"random phase approximation" (RPA). Then, taking into account Eq. (JTUJ) one obtains [17,18] 

U RPA {q,z) =Il < - \q : z) + n T (q,z). (28) 

n (0) (^) = ^%- (29) 



*r, ,_ / d 3 p /S(p~q/2)-/S(P + q/2) 



and 

n ' (9 ' Z) = J teTe(p^-k/2)-e(p'+k/2) (30) 

Taking into account the above consideration, we now modify the RPA approximation on the generalized (MRPA) 
approximation, by introducing the change in H T (q,z) lj30j) operator. The function ff d (fTT]) and the spectrum of the 
one-particle excitations e(p) for the ideal Bose gas are changing for the function f T (p) and the energy E*(p) lfl4|) in 
the self-consistent Hartree-Fock approximation. As this takes place, we observed that due to the ruptured character 
of the one-particle spectrum at the point p = 0, the function H^(q, z) preserves its earlier form lj29|) . Then 

U MRPA (q,z) = n^(q,z) + n T MRPA (q,z), (31) 
n T , s f d Z V / T ( P -q/2)-/ T (p + q/2) 

lL MRPA [q, z)-J ^3 hz + E * {p _ k/2) _ E « {p + k/2 y W 

Considering further the case of a low temperature in ((19)1 . one can omit the part n^ fi j PA (g, z) ([32)1 in U MRPA (q 7 z). 
In this case, from l[25 |l -p7 | . (|29| it directly follows that for the temperatures T <C A 

where the spectrum of the collective excitations is determined by the equality 

hu;( q ) = {(e(q)) 2 + 2e(q)nu(q)} 1/2 . (34) 
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The relation (|34|) completely coincides with the spectrum of the "bogolons" for the interaction potential, where its 
dependence on the wave vector [26,27] provides the fulfilment of the Landau superfluidity condition. Since, in the 
temperature interval under consideration, the density of the particles in the "codensate" n is close to the total 
density n, we can change no to n in f34|) . Then, inserting f33|) in f22|) and taking into account the determination of 
the structure factor [1,23], we find [17,18] 

nS {q ) = ^cth\^m (35) 



fuv(q) { 2T J 

The above equation is the generalization of the Feinman formula [28] for the connection between the static structure 
factor and the spectrum 

M<?) = S(q)/S(q), (36) 

which is valid for the case hw(q) ^> T. 

In the opposite case of huj(q) <C T (and no ~ n, e.g., T <C To), we have 

h' L LU' L {q) q^Q nu(0) 

Relation (|37|l corresponds to the general result for the systems with the short-range interaction potential (1) [1] 



1 fdV\ 

]imS(q)=nTK T , K T = --[ — ) . (38) 



w J ' V \dP J 

Here Kt is the isothermal compressibility of the system. Comparing |37|) and (|38| and taking into account that the 
spectrum of the "bogolons" in the region of the small wave vectors has the form 



h^ q ) Sj q . Sj = I'^lV 2 (39) 



we conclude that the value St characterizes the isothermal sound velocity. In addition, on the basis of Eq. (|35|) for the 
static structure factor we can determine the free energy F of weakly non-ideal Bose gas for the temperature T -C A 
by using the general relation [1] 



F = F >« + l u (0)nN - \n ^ u(q) + \n ^ J 



u(q)S x (q)d\, (40) 



F id = T^2ln[l-exp(-e(q)/T)}. (41) 

9 

Here F ld is the free energy of the ideal Bose gas for T < Tq [26] and 5a (<?) is the static structure factor for the system 
with the interaction potential Xu(q). Inserting (|35|l in (jJOj we find [17,18] 

F = \u{0)nN - ~ ^{£(9) + nu(q) - {e(qf + 2e(q)nu(q)) 1/2 } + {1 - exp[-^(g)/T]} . (42) 

9 9 

Relation f42|) completely corresponds to the results obtained in [5,24,27]. 



IV. CONCLUSIONS 



In the present work we developed the theory of weakly non-ideal Bose gas, below the condensation temperature, on 
the basis of the self-consistent Hartree-Fock approximation. It was found that the spectra of the collective excitations 
and the one-particle excitations are distinct, which is contrary to the theory based on the C-number representation 
for the operators and a . It is shown that one-particle branch of excitations has a gap. Both spectra satisfy the 
Landau criterion of superfluidity. It must be noted that in the early works of Landau, and in the works of Bogolyubov, 
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they considered the possibility of the existence of the gap in the spectrum of Bose system for T < Tq. However, they 
omitted this suggestion because the phonon branch was absent. 

It follows now from the present paper, that the collective phonon-roton spectrum and the one-particle spectrum with 
a gap actually coexist. On the basis of the calculation of the one-particle distribution function and "density-density" 
Green function, in the framework of the self-consistent Hartree-Fock approximation for the weakly non-ideal Bose 
gas, we can establish the following: 

1) This system possesses two branches of excitations - the one-particle branch and the collective branch, each of 
them satisfying the Landau condition of superfluidity. 

2) In the region of small wave numbers, there is the gap in the spectrum of the one-particle excitations which is 
conditioned by the presence of the "condensate". 

3) The spectrum of the "bogolons" corresponds to the phonon-roton excitations that are observed in the experiments 
on the neutron nonelastic collisions [29,30]. 

Therefore, even a weak inter-particle interaction leads to the drastic differences in the description of the one-particle 
distribution function and of the excitations of the "overcondensed" particles. We have emphasizes that (as in the case 
of C-number representation of the operators Oq and ao) the application of the Hartree-Fock approximation to the 
calculation of the Green functions for Bose gas in the temperature region T < Tq cannot be absolutely rigorously 
theoretically justified and there is a need here for further experimental examination. The principal difference between 
the results of this paper and the traditional approach based on the C-number representation for aj and ao, is the 
appearance of the gap in the spectrum of the one-particle excitations. It follows from the results of this work that 
the gap cannot manifest itself in the "density-density" Green function (at least in the Hartree-Fock approximation) 
and therefore it cannot be seen in the experiments on neutron scattering in the superfluid Helium [29,30]. However, 
such possibility cannot be excluded in the experiments on Raman light scattering. Moreover, in [31], where such 
experiments are described, there is the direct indication of the existence of the gap. 



Acknowledgment 

The authors thank Yu. A. Kuharenko for the useful discussion. V.B. and S.T. express gratitude to the Netherlands 
Organization for Scientific Research (NWO) for support of their investigations on the problems of statistical physics. 



[I] R. Balescu, Equilibrium and nonequilibrium statistical mechanics, v.l, J. Wiley and Sons, Inc. (1975) 
[2] L.D. Landau, Zh. Exp. Teor. Fiz. 11 (1941) 592 

[3] L.D. Landau, Phys. Rev. 60 (1941) 356; J.Phys.USSR 5 (1941) 71 
[4] L.D. Landau, J.Phys.USSR 11 (1947) 91 

[5] A. A. Abrikosov, L.P. Gorkov, I.E. Dzyaloshinskii, Quantum Field Theoretical Methods in Statistical Physics, 2ed., 
Pergamon, (1965) 

[6] N.N. Bogolubov, J.Phys.USSR 11 (1947) 23 

[7] N.N. Bogolubov, Lectures on quantum statistics, Kyiv "Radianska Shkola"1949; Selected works v. 2 (1969). 
[8] S.T. Belyaev, Zh. Exp. Teor. Fiz. 34 (1958) 417, 433. 
[9] N.M.Hugenholtz and D. Pines, Phys.Rev. 116 (1959) 489 
[10] W.H.Bassichis and L.L.Foldy, Phys.Rev. 133A (1964) 435 

[II] H.Stolz, Physica A86 (1977) 11 

[12] I.C. Landman, S.E. Peschanov, FR. Ulinich, Preprint of the Institute of Atomic Energy N248 97/1 (1989) 
[13] I.M. Yurin, E-print archives, quant-ph./0310115 

[14] I.M. Yurin and S.A. Trigger, E-print archives, cond-mat.supr-con./0906.0755 
[15] S.A. Trigger and P.P.J. M.Schram, Physica B 228 (1996) 107 

[16] S.A. Trigger, I.M. Yurin, Extensions of the Bogolyubov's theory of weakly non-ideal bose gas, Modern problems of 
theoretical and mathematical physics, Bogolyubov Kiev conference, Book of abstracts, 36 (2009) 
[17] V.B. Bobrov and S.A. Trigger, Physica A 170 (1990) 261 

[18] V.B. Bobrov, S.A. Trigger, Yu.P. Vlasov, Zh. Eksp. Teor. Fiz. 102 (1992) 107 (Sov.Phys. JETP 80 (1995) 853) 

[19] V.B. Bobrov, S.A. Trigger, P.P.J. M. Schram, Zh. Eksp. Teor. Fiz. 107 (1995) 1527, [Sov.Phys. JETP 80 (1995) 853] 

[20] D.A. Kirznitz, The field method of the many-body theory. Moscow, Fizmatgiz (1963) 

[21] L.P. Kadanoff, G. Baym, Quantum statistical mechanics, W.A. Benjamin, Inc. N.-Y. (1962) 

[22] Yu.L. Klimontovich, V.P. Silin, JETP 23 (1952) 151 

[23] N.H. March, N. and M. Parinello, Collective Effects in Solids and Liquids, Bristol: Hilger, 1982 

[24] E.M. Lifshitz, L.P. Pitaevsky, Statistical Physics, part 2 Thoery of condensed matter, Nauka, Moscow (1978) [in Russian] 
[25] W.-D. Kraeft, D. Kremp, W. Ebeling, G. Ropke, Quantum statistics of charged particle systems. Berlin: Akademie- 
Verlag (1986) 



8 



[26] Yu. A. Zerkovnikov, Theor. Math. Phys. 4 (1970) 119 [in Russian] 

[27] A. I. Akhiezer and S. V. Peletminski, Methods of Statistical Physics, Pergamon, New York, (1981) 
[28] R.Feynman, Phys.Rev. 94 (1954) 262 

[29] R.A.Cowley and A.D.B.Woods, Can.J.Phys. 49 (1971) 177 

[30] N.M. Blagoveshchensky, E.B. Dokukin, J. A. Kozlov, V.A. Parfenov, Pis'ma v ZhETF 31 (1980) 7 
[31] T.J.Greytak, R. Woerner, J.Yan, R.Benjamin, Phys. Rev. Lett. 25 (1970) 1547 



